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to : ABSTRACT 

Stationary solutions to the problem of particle acceleration at shock waves in 
s ! ' the non-linear regime, when the dynamical reaction of the accelerated parti- 

Q\ ' cles on the shock cannot be neglected, are known to show a prominent energy 

IT) '. 

pvi , flux escaping from the shock towards upstream infinity. On physical grounds, 

^. . the escape of particles from the upstream region of a shock has to be expected 

f**^ ■ in all those situations in which the maximum momentum of accelerated par- 



ticles, p m ax, decreases with time, as is the case for the Sedov- Taylor phase 

of expansion of a shell Supernova Remnant, when both the shock velocity 
•i— I , 
^^ . and the cosmic ray induced magnetization decrease. In this situation, at each 

^j ■ 

C^ ■ time t, particles with momenta larger than p ma x(t) leave the system from 

upstream, carrying away a large fraction of the energy if the shock is strongly 
modified by the presence of cosmic rays. This phenomenon is of crucial im- 
portance for explaining the cosmic ray spectrum detected at Earth. In this 
paper we discuss how this escape flux appears in the different approaches to 
non-linear diffusive shock acceleration, and especially in the quasi-stationary 
semi-analytical kinetic ones. We apply our calculations to the Sedov- Taylor 
phase of a typical supernova remnant, including in a self-consistent way par- 
ticle acceleration, magnetic field amplification and the dynamical reaction on 
the shock structure of both particles and fields. Within this framework we 
calculate the temporal evolution of the maximum energy reached by the ac- 
celerated particles and of the escape flux towards upstream infinity. The latter 
quantity is directly related to the cosmic ray spectrum detected at Earth. 

Key words: acceleration of particles - shock waves 
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1 INTRODUCTION 

Supernova Remnants (SNRs) have long been suspected to be the main sources of galactic 
cosmic rays with energies up to the so called "knee" (£~3x 10 15 eV). The existence of a 
flux of energetic particles escaping the accelerator from the upstream region is essential if 
this paradigm is indeed realized in Nature. The commonly accepted view of a typical SNR 
evolution is as follows. During the initial stage after the supernova explosion, the ejecta are in 
free-expansion with a high velocity ~ 10 4 km/s. Non-linear acceleration is expected to be at 
work: efficient particle acceleration is associated with shock modification and magnetic field 
amplification in a complex chain of causality. During this stage, the maximum momentum 
achieved by the a ccelerated particles increases with time up to p ma x ~ 10 16 eV/c (see e.g. 



Blasi ct al. 



(J2007I ) and references therein). When the mass of the swept-up material becomes 



comparable with the mass of the ejecta, the Sedov- Taylor phase begin s: the shock slows down 



and m agnetic field damping is expected to become more effective ( iPtuskin Sz Zirakashvili 



(120051 )). As a consequence, the maximum momentum of particles that can return to the 
shock is expected to decrease. At each time t, particles with momentum exceeding the 
current p m ax(t) do not make it back to the shock and leave the system from upstream: this 
is the so called "escape flux" . 

The escape flux plays an essential role in the formation of the cosmic ray spectrum de- 
tected at Earth. If there were no escape from upstream during the Sedov- Taylor phase, all 
particles accelerated in a SNR would be advected downstream and undergo severe adiabatic 
energy losses before being injected into the ISM: in this case S NRs could not b e resp on- 
sible for accelera ting cosmic rays with energie s around the knee ( 1 Caprioli et al.l ( 120081 ). in 
preparation, and lPtuskin fc Zirakashvili! (120051 )). 

While being of crucial importance during the Sedov- Taylor phase, it is clear that escape 
cannot occur during the free expansion phase of a SNR: the only particles that can leave the 
system from upstream are the ones at the highest energies, but, since p m ax is now increasing 
with time, no appreciable energy flux can escape the accel erator during this phase. 



Kin etic approaches to non-linear particle acc eleration (jMalkovl (119971 ): 



(120011 ) : iBlasil (12002 



Malkov fc Drury 



20041 ) : lAmato fc Blasil (120051 )) allow us to calculate the spectrum and the 



spatial distribution (including the absolute normalization) of the particles accelerated at the 
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shock, even in the case when the diffusion coefficient is the result of magnetic field amplifica- 



tion b y streaming instability induced by the accelerated particles themselves ( jAmato fc Blasi 



(200 6f)). These appr o aches , as well as others that have appeared in the literature (for instance 



Berezhko fc Ellison! (119991 )). are based on the assumption that the acceleration process (and 
the shock modification) may be assumed to reach some sort of stationarity. In all these cases 
the calculations show that the shock is strongly modified by the presence of cosmic rays, 
and that the spectra are concave, with a slope at momenta close to the maximum momen- 
tum Pmax which is flatter than p~ A . All these calculations, independently of the techniques 
used to solve the equations, predict an escape flux of particles (and energy) towards up- 
stream infinity: the shock becomes radiative, which is one of the very reasons why the shock 
modification becomes effective (namely the total compression factor becomes larger than 4). 

It is worth recalling that t he assumption o f stati onarity was widely used also in the 
framework of two- fluid models (jDrury fc Volki (11981a! Jbl)). and its relevance in the case of 
mo dified shocks, le a ding to particle spectra harder than p~ A at high energies, was discussed 
by iMalkov fc Volki (119961 ) in the context of the so-called renormalized two-fluid models. It 
should be stressed, indeed, that the escape flux can arise both in the linear and non-linear 
regime of particle acceleration. The important difference between the two regimes, however, 
is that within the former the escaping particles carry away a negligible energy flux, while 
for modified shocks the bulk of the energy in the form of cosmic rays is stored in the highest 
energy particles and leaves the system with them. 

The main weak point of the assumption of quasi- stationarity is that it leads to artificial 
escape fluxes in situations, such as the free expansion phase of a SNR, where escape would 
be unphysical (at least in the idealized cases of plane and perfectly spherical shocks). In 
this case, the appearance of an escape flux signals for the need for a fully time dependent 
approach. In other words, what appears as an escape flux is the energy which is channeled 
into particles with ever increasing momentum. During the Sedov- Taylor phase, as emphasized 
earlier, escape is in fact possible if the magnetic field in the shock region decreases. In 
this case, the leakage of particles towards upstream allows the system to relax to a quasi- 
stationary situation. 

A number of ti me-dependent stud i es of n on- linear shock dynamics have been presented in 



Berezhko fc Volk 



the lit e rature (e.g. iFalle fc Giddingd (J19871 )). in particular for SNRs (e.g. 

( 19971 ); iKang fc Joned ( 20061 ) and references therein). All these calculations aim at solving 

the transport equations and the equation of fluid dynamics numerically. Most of these papers 
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have however concentrated on the spectrum of accelerated particles rather than on the 
escape flux. In fact, in many of these papers the magnetic field is introduced by hand and 
is not a function of time, thereby not leading to a decrease with time of the maximum 
momentum of accelerated particles. In other cases the boundary conditions are imposed at 
upstream infinity, therefore the escape flux, even if it is present, would appear in the form 
of concentration of high energy particles at large distances from the shock. In these cases a 
recipe is needed to label these particles as escaping the system (see for instance ?). 



The identification o 



the escape flux is all but a mathematical detail: as discussed by 



Ptuskin fc Zirakashvilil (120051 ) . the escape flux is crucial for establishing a connection between 
SNRs and the origin of cosmic rays, in that it allows the escaping particles to avoid the 
adiabatic losses related to the expansion of the supernova shell, and therefore become cosmic 
rays at the knee. In the absence of escape towards upstream infinity at the beginning of the 
Sedov- Taylor phase, SNRs cannot be the sources of cosmic rays up to the knee. 

The c ontribution of the escape flux to the galactic cosmic ray spectrum was first esti- 



mated by lPtuskin fc Zirakashvilil (120051 ). These authors included the magnetic field damping 
and the consequent decrease of p ma x during the Sedov- Taylor expansion, but in their study 
the shock structure and particle spectrum are fixed rather than self-consistently calculated. 

In this paper we present the first attempt at carrying out a calculation of the escape flux 
within the framework of a kinetic approach to non-linear shock acceleration, self-consistently 
including particle acceleration, magnetic field amplification and the dynamical reaction of 
both on the shock. We focus our attention on shocks in shell type SNRs, discussing the 
physical meaning of the escape flux during the different phases of their evolution. In par- 
ticular the flux of energetic particles leaving the remnant during the Sedov- Taylor phase is 
explicitely computed and its phenomenological implications for the origin of cosmic rays are 
discussed. 

The paper is organized as follows: in $2] we discuss the implications of the assumption 
of stationarity of the acceleration process. In $3] we discuss the escape flux based on the 
most general version of the conservation equations. In §H we discuss how the escape flux 
is connected to the existence of a maximum momentum in the distribution of accelerated 
particles. In £J5] we apply our calculations to the different stages of evolution of a supernova 
remnant. We conclude in $6l 
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2 THE ASSUMPTION OF STATIONARITY 



The standard solution of the stationary transport equation 



u\x) 



d f(x,p) 
dx 



d_ 

dx 



D{p) 



df(x,p) 
dx 



pdudf 
3 dx dp 



(1) 



leads, in the test-particle regime, to the well known power-law spectrum of accelerated 
particles f(p) oc p~ a , with a = 3r/(r — 1) where r is the compression factor at the shock. 

The power law extends to infinitely large momenta. Since for ordinary non-relativistic 
gaseous shocks r < 4 (namely a > 4), the total energy in the form of accelerated particles 
remains finite. This solution is found by imposing as boundary condition at upstream infinity 
(x = — oo) that /(— oo ) = and d f {—oo) / dx = 0. 

If the boundary condition f(x = x ) = is used, instead, at some finite distance x < 
upstream, the solution of the transport equation is easily calculated to be 



f{x,p) 



Up) 



exp 



U\X 



where 

f (p) = Kexp 



l-exp(gg) I \D(p), 
3u\ f'P dp' 



— exp 



MP), 



U\ 



i - «a Jpin) V' 1 - exp (») 
In the case of Bohm diffusion, D(p) = Do(p/m p c) and one obtains: 

, ■ , , , t 3m i f p d P' 1 

Jo(P) = K exp 



(2) 



(3) 



(4) 



(/, 11-2 J Pinj p' l_ exp ^_£i^ 

xo\uim p c/ Dq. Now one can show that for p <^ p*, fo(p) oc {p/p*)~' Zr ^ r ~ l \ with 



3r p 

r— 1 p» 



The quantity 



where p* 

r = U1/112, the standard result. However, for p ^> p*, fo(p) oc exp 

Pmax = P*{f — l)/3r plays the role of maximum momentum of the accelerated particles. 

This simple example shows how a maximum momentum can be obtained in a stationary 
approach only by imposing the boundary condition at a finite boundary. Physically this 
corresponds to particles' escape, as shown by the fact that the flux of particles at x = x$ is 

df(x ) uifoip) fuix s 



<p(x ,p) = Ui/(a? ,p) - D(p)- 



exp 



< 0. 



(5) 



dx i_ exp (W) \D(p) / 

The fact that 4>(xq,p) < shows that the flux of particles is directed towards upstream 
infinity. Moreover, the escape flux as a function of momentum, <ft(xo,p), is negligible for all 
p with the exception of a narrow region around p m ax'- only particles with momentum close 
to Pmax can escape the system towards upstream infinity. The escape flux as a function of 
momentum is plotted in Fig. [1] for two values of the compression factor, r = 4 (solid line) 
and r = 7 (dashed line). The normalizations are arbitrary, since the calculations are carried 
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Figure 1. We plot the escape flux (f>(xo,p) as a function of momentum. The curves refer to two different values of the shock 
compression ratio: r = 4 (solid line) and r = 7 (dashed line). The computation is carried out in the test-particle regime. The 
x-axis is in units of the reference momentum p„ = r/(r — l)p m ax, while units along the y-axis are arbitrary 



out in the context of test particle theory. The latter value of r cannot be realized at purely 
gaseous shocks, but we have adopted this value to mimic the effect of shock modification, 
which leads to total compression factors larger than 4. 

The escape phenomenon is basically irrelevant in the test-particle regime, because of the 
negligible fraction of energy carried by particles with p ~ p m ax, but it becomes extremely 
important in the calculation of the shock modification induced by accelerated particles. For 
strongly modified shocks, the slope of the spectrum at high energies is flatter than p~ A and 
the fraction of energy that leaves the system towards upstream infinity may dominate the 
energy budget. This is the escape flux which appears in all approaches to cosmic ray modified 
shocks. 

In the context of kinetic calculations of the shock modification in the stationary regime, 
the escape flux appears however not as a consequence of imposing a boundary condition at 
a finite dista nce upstream, but r a ther as an apparent violation of the equation of energy 
conservation (IBerezhko fc Ellison! (119991 )). that requires the introduction of an escape term 
at upstream infinity. In the next section we discuss this effect, which reveals the true nature 
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of the escape flux, as related to the form of the conservation equations and the assumption 
of stationarity. 



3 CONSERVATION EQUATIONS AND ESCAPE FLUX 

In this section we rederive the conservation equations for cosmic ray modified shocks in their 
general form, in order to emphasize the mathematical origin of the escape flux. 

The time dependent conservation equations in the presence of accelerated particles at a 
shock can be written in the following form: 
dp d(pu) 



d_ 
dt 



dt 

d(pu) 

dt 

Pa 



- 2 pu 



7,-lJ 



dx 

d_ 

dx 

d_ 

dx 



pu 



P„ + Pr + P 



w 



2 7,-lJ 



d 



u— [P c + P w ] + TE W 



(6) 

(7) 
(8) 



Here P g , P c and Pyy are respectively the gas pressure, the cosmic ray pressure and the 
pressure in the form of waves. E\y is the energy density in the form of waves and V is the 
rate at which the background plasma is heated due to the damping of waves onto the plasma. 
The rate of change of the gas temperature is related to TE\y through: 



dP„ dP„ _, du . lXinrn 

The cosmic ray pressure can be calculated from the transport equation: 



df(t,x,p) df(t,x,p) 

+ u\x) 



dt 



dx 



d_ 

dx 



D(x,p) 



df(t,x,p) 
dx 



pdf(t,x,p) du(x) 
3 



dp 



</.;■ 



(9) 



(10) 



where we put u(x) = u(x) — vw{x) and v\y{x) is the wave velocity. For our purposes here 
we are neglecting the injection term. 

Multiplying this equation by the kinetic energy T(p) = m p c 2 [^i — l), where 7 is the 
Lorentz factor of a particle with momentum p, and integrating the transport equation in 
momentum, one has: 

du 



dE c d{uE c 



dt 
where 

E r 



dx 



d_ 

dx 



D 



.dEc 
dx 



P 



dx ) 



f°° f 00 4:71 

/ dp4.iip 2 T(p)f(p) and P c = / dp—pv(p)f(p) 

Jo Jo o 



(11) 



(12) 



are the energy density and pressure in the form of accelerated particles. Moreover we intro- 
duced the mean diffusion coefficient: 
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D(x) 



£> AmpT{p)U 



dx 



(13) 



The only assumption that we made here is that f(p) — > for p — > oo. We shall comment 
later (§3| on what would happen if the spectrum were truncated at some fixed p max , instead. 

Introducing the adiabatic index for cosmic rays 7 C as E c 
Eq. E] as 



— - H 

<9t cte 



7c^-P c 
7c -1 



9a; 



dx 



+ u- 



dPc 
dx 



Pc/(lc — 1); w e can rewrite 

(14) 



and use it to derive udP c /dx. In this way Eq. [H] becomes: 

d_ 
dt 



1 2^ Pg 
-pu + y - 

2 lg~ 



+ E C 



d_ 
dx 



-pu° + 

2 lg 



IgPgU IcPcU 



+ 



7c 



+ 



d_ 

dx 



D{X) ^X~ 



-%■ 



d_Pc 
dx 



dPi 



— it- 



ir 



dx 



+ TE 



ii'- 



(15) 



At this point we can make use of the equation describing the evolution of the waves: 



dEw dF\Y 



dPw rp r rp 

u— h crh w — 1 L w , 



dt ' dx dx 

where a is the growth rate of waves, integrated over wavenumber and Fy/ is the flux of energy 
in the shape of waves. These quantities can be calculated once it is known how particles 
with given momentum p interact with waves with wavenumber k and how the turbulence 
energy is shared between the kinetic and the properly magnetic contributions. Substituting 

into Eq. [15] we get: 
d_ 

1ft 



1 2^ Pg 
-pu + y — 

2 7 ff -l 



+ E C +E: 



It' 



d_ 
dx 



2 7 S 



IgPgU IcPcU 



+ 



7c 



+ *1 



11' 



d_ 
dx 



D(x) 



dEc 
dx 



%■ 



dPc 

dx 



+ a3 



w ■ 



(17) 



In the case of Alfven waves resonan t with the Lar mor radius of the accelerated particles, 
one has v w = v A = Bji^Txpfl 2 and Jskillinj (119751 1): 



oE 



IV 



VA 



d_Pc 

dx 



(18) 



so that the energy conservation equation reads 



d_ 
dt 



-pu 2 4 



lg 



+ E C + E< 



w 



d_ 
dx 



1 ■} IqPqU IcPcU „ 

2 Jg - 1 7c - 1 



+ 



d_ 

dx 



D(x 



dE c 
dx 



(19) 



In the general case of waves other than resonant Alfven waves, Eq. [18] does not hold 
and even the equality vw = va may be questionable: as a consequence one cannot use the 
standard Eq. [19], while Eq. [T7] is still correct. However, in order to be able to solve the 
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problem, the detailed form of the cosmic ray transport equation, as well as an expression 
analogous to Eq. [181 relating the growth of the wave energy to the cosmic ray dynamics, 
are still needed. It is worth stressing, in fact, that the magnetic turbulence has often been 
proposed to not show up in the fo rm of stan dard resonant Alfven waves, but rather as 



non-resonant purely growin g modes 



phenomenological model of 



Bell 



Bell fc Lucek 



120041). or as generic "magnetic structures" in the 



(120011 ) . This also causes the connection between 
F w and Pw to be generally different from the standard F w « 3uPw, thus the contribution of 
these waves to the energy conservation equa t ion d oes not necessarily lead to Eq. [19], which 



was neverthe 
suggested by 



ess used by 



Vladimirov et al. 



(120061 ) in the implementation of the scenario 
Bell fc Lucekl (120011 ). 

In the following we limit ourselves to the case of Alfven waves, which interact resonantly 
with particles, since in this case the calculations are all well defined. 



Notice that in the stationary regime, Eq. [HI integrated around the subshock leads to 



7c 



-uP r - D 



dE, 



constant, (20) 

7 C — 1 dx 

because of the continuity of the cosmic ray distribution function. On the other hand, Eq. [TH] 
(again in the stationary case) , once integrated around the shock, leads to conclude that: 

-pu 3 + 7g gU +E W = constant. (21) 

2 7g~l 

In other words, at the subshock the energy fluxes of the gaseous and cosmic ray components 
are conserved separately. This is what is usually meant when we refer to the subshock as an 
ordinary gas shock. In the following we use the stationary version of Eq. 



d_ 

dx 



IgPgU IcPcU 



-pu* + -^- + 

2 Jg ~ 1 lc 



1 



+ F, 



IV 



D(x 



dE c 
dx 



0. 



(22) 



4 ESCAPE FLUX AND THE NEED FOR A P MAX 

In non-linear theories of particle acceleration the need for a maximum momentum is dic- 
tated by the fact that the spectrum at large momenta becomes harder than p~ 4 , so that 
in the absence of a high p cutoff the energy content of the accelerated particle distribu- 
tion would diverge. Before this happens the dynamical reaction of the accelerated particles 



would i nhibit further accele r ation. In most ap p roaches to non-linea r calcu lations ( jMalkov 



(1997) 



Malkov et al. 



(12Q00|); 



Blasi (2002 



20041 ); Berezhko k Ellison! (119991 )). the maximum 



momentum is a given parameter, taken together with the assumption of stationarity of the 
acceleration process. The transport equation is then solved between the shock and upstream 
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infinity. Both in the downstream region and at upstream infinity one has Ddf/dx 



0. 



Moreover, at upstream infinity there are no accelerated particles (P c = 0) so that Eq. 1221 
becomes: 



2 la -1 



lgPg,2U2 . 1cP C U2 



+ 



+ F } 



w 



1 3 , 1g P g,o u o 
2 P ° U ° + ~^T 



, .. .„..., , (23) 

None of the calculations of particle acceleration at modified shocks carried out so far satisfies 
Eq. [23] unless it is completed with an escape flux F esc such that 



\" 2Ui+ 7,-1 



lgPg,2U2 , 7c-Pc«2 



+ 



+ ^1 



W 



~Pou Q + 

2 lg 



IgPgflUo 



-F P . 



7c -l " -I 1 -"" -„-l '- (24) 

Unfortunately, as showed above, this apparently harmless step is inconsistent with f(p,x) 
being a solution of the time-independent transport equation. In fact this is not surprising 
since, as we stressed above, the solution of such equation cannot be characterized by a finite 
Pmax when the boundary condition of vanishing f(p,x) and df/dx is imposed at upstream 
infinity. 

On the other hand, as also discussed by iMalkov fc Volkl (119961 ) in the context of the so- 
called renormalized two-fluid models, the requirement that there is a maximum momentum 
leads to the appearance of an additional term on the right hand side of Eq. [HI which is 

+ 7T-T 47ip 3 T(p)f( P ,x) 

Sax L J p=o 



(25) 



This reflects in an additional term in Eq 
r° + 1 du 



/O.X z t7rp max l [Pmax ) J [Pmax j X ) - 

-oo 6 dX 

^P 3 max T (Pmax)f(Pmax, x) + («2 ~ «l) 



(26) 



0_ A ldU A 3 

-oo o ax 



Air 



x ■*■ (Pmax ) JO \Pmax ) ■ 



Since du/dx < in the precursor, this term is negative and numerically coincides with the 
escape flux F esc . 

As an alternative to obtaining the escape flux in this way, one could derive it by assuming 
that there is no p max imposed by hand, and that a maximum momentum results from 
imposing th e boundary condi t ion at a finite distance upstream, rather than at upstream 



infinity (e.g. 



Vladimirov et al. 



(200 



3))- 



In this second case Ddf/dx does not vanish at the upstream boundary and an escape 
flux appears in a natural way, namely 



<p eS c = Uof{xo,p) - D 



dx 



-D 



J X=XQ 



dx 



<0, 



(27) 



X=XQ 



and the energy escape flux F esc is related to 4> esc through 
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in a ,'('■ 



Fesc = 4vr p 2 dp <f) esc {p) T{p) . (28) 

J Pinj 

In other words, no artificial escape flux needs to be introduced if a boundary condition 
is imposed at a finite distance upstream, or if, as an alternative, the fully time dependent 
solution of the problem can be found. 

In next section we explore the consequences of the existence of an escape flux for the 
origin of cosmic rays in supernova remnants. 

5 PHYSICAL MEANING OF THE ESCAPE FLUX FOR SUPERNOVA 

REMNANTS 

The acceleration process in supernova remnants is expected to work in qualitatively different 
ways during the free expansion and the Sedov- Taylor phases. Here we restrict our attention 
to the propagation in a spatially uniform interstellar medium. During the free expansion 
phase the velocity of the shell remains constant and the maximum momentum grows in 
time in a way that depends on the growth of the turbulent magnetic field in the upstream 
region. During this phase particles cannot escape. Nevertheless the standard approaches to 
the calculation of the shock modification would lead to predict an escape flux, a symptom of 
the need to carry out fully time dependent calculations to treat this expansion regime. The 
lack of particles' escape implies an increase in the maximum momentum of the accelerated 
particles. This trend ends at the beginning of the Sedov- Taylor phase, when the inertia of 
the swept up material slows down the expanding shell. Physically, this is the reason why 
we expect that the highest energies for particles accelerated in SNRs are reached at the 
beginning of the Sedov- Taylor phase. 

During this phase, the shock velocity decreases and the magnetic field amplification 
upstream, as due to streaming instability, becomes less efficient. The generat ion of magn etic 



turbulence via st reaming in stabilities may proceed through either resonant (jBelll (Il978al jbl)) 



or non-resonant (JBelll (12004 ) ) coupling between particles and waves and the t wo channels are 



l ikely t o dominate at di f ferent times in the history of the supernova remnant (IPelletier et al.l 



(J20061 ): lAmato fc Blasil (120081 ) ). as discussed below. 



This general picture leads to a maximum momentum that decreases with time and to 
particles' escape towards upstream infinity: particles of momentum p max {ti) do not make it 
back to the shock at a time t% >t\. In other words, during the time interval between t\ and 
t 2 , particles with momentum between p max (ti) and p max (t 2 ) escape from the system. This 
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happens at any time, and a net flux of particles (and energy) towards upstream infinity is 
realized. At any given time t the spectrum of particles that escape is highly peaked around 
Vmaxif) (see Fig. Q] for the test-particle case). The spectrum of accelerated particles that 
is confined in the accelerator and advected towards downstream is cut off at a gradually 
lower maximum momentum, and this should reflect in the spectrum of secondary radiation, 
especially gamma-rays. The particles trapped downstream will also eventually escape the 
system, but at the time this happens they will have been affected by adiabatic losses due 
to the expansion of the shell. Therefore this part of the escaping flux, which will reflect the 
history of the remnant, is suspected to play a particularly important role only at the lowest 
energies in the cosmic ray spectrum at earth. 

The flux of high energy cosmic rays, close to the knee region, as we see below, is mainly 
generated during the early Sedov- Taylor phase and is made of particles escaping the accel- 
erators from upstream. The actual flux of diffuse cosmic rays observed at the Earth results 
from the integration over time of all the instantaneous spectra of escaping particles, each 
peaked at p max (t), and from the superposition of the supernova explosions that could con- 
tribute. This integration is affected by the accelerator properties, by the dynamics of the 
expanding shell and by the damping processes that ma y affect the way the magnetic field is 
amplified by streaming instabilities at any given time (IPtuskin fc Zirakashvilil (120051 )). 

There is also another implication of the line of thought illustrated above: the spectrum of 
particles that escape, as integrated over time during the Sedov- Taylor phase, does not need to 
be identical to the spectrum of particles advected towards downstream. However, the latter 
are the particles which are responsible for the production of secondary radiation (radio, 
X-rays, gamma rays): the concave spectra predicted by the non-linear theory of particle 
acceleration and to some extent required to explain observations, might not be reflected in 
a concavity of the spectrum of escaping particles. 

An estimate of the scalings of the relevant quantities during the ST phase can be found 
as follows. The radius and velocity of the expanding shell can be written as: 

R sh (t) = 2.7 x 10 19 cm (^V \%l (29) 

\ TIq / 

V sh (t) = 4.7 x 10 8 cm/s (^i) t$ 5 , (30) 

where E 5 i is the kinetic energy of the shell in the free expansion phase in units of 10 51 erg 
and no is the number density of the plasma upstream in units of cm -3 . Here we assumed the 



E max (t) = 3.8 x 10 4 5B, G (t) ( ^1 ) i"V5 GeV (31) 



/ FT \ V 10 

*fl(t) = 65 n 1 / 4 < G (^1) t-/ 10 ^) 1/2 ^ . (32) 
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standard Sedov- Taylor time-scaling of R s h and V s h, but the reader should bear in mind that 
the adiabatic solution may be affected by the fact that in this phase the shock is radiating 
energy in the form of cosmic rays. The maximum energy is estimated by requiring that the 
diffusion length upstream equals some fraction (say 10%) of R s h(t). If the diffusion coefficient 
is assumed to be Bohm-like and the magnetic field close to the shock is SB(t), one obtains: 

'£5l\ 2/5 

n 
The magnetic field in the shock vicinity is amplified by streaming instability, induced by the 

accelerated particles both resonantly and non-resonant ly. Let us introduce the acceleration 

efficiency as a function of time: £ c (£) = P c (t) / (poVshit) 2 ) . In terms of £ c , the strength of 

the resonantly amplified magnetic field at the saturation level can be estimated as: 5B 2 = 

8ttpqV 2 ^ c /Ma (Ma is the Alfven Mach number), which leads to: 

- rr \ 1 / 10 
n 
In a similar way, the strength of the field in the case of non-resonant amplification can be 

estimated from 5B 2 = 2tt po{V s h(t) 3 / c)£ c (i) and leads to: 

' rr \ 3 / 10 

■&51N —fi 

n 

In general the two channels of magnetic field amplification work together but the non- 
resonant channel dominates at earlier times and leads to stronger magnetic field amplifica- 
tion. 

The maximum momentum in the two cases is as follows: 

■ rr \ V2 

n 
in the resonant case, and 

/rr \ 7 /io 
E max (t) = 7.3 x 10 6 (£^J nT Ut) 1/2 C A0 GeV ( 35 ) 

in the non- resonant regime. 

In the naive assumption that the acceleration efficiency is constant in time, we see that 
E ma x(t) scales with time as £~ n / 10 at earlier times and as £ -1 / 2 at later times, when resonant 
scattering dominates. In actuality the scalings will be more complex because of the non-linear 
effects (especially the formation of a precursor upstream) induced by accelerated particles, 
which also lead to a time dependence of £ c (0- 

As discussed in the previous sections, it is not clear how to describe the non-resonant 
waves in the context of the conservation equations. A calculation of the dynamical effect of 



SB(t) = 198 n 1 / 2 ( ^1 ) ^fo ^1/2 ^ (33) 



E max (t) = 2.5 x 10 6 (^i) V2 uTbZg Ut) 1/2 tlf GeV, (34) 
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these modes on the shock is therefore not reliable at the present time. For this reason, here 
we confine ourselves to the investigation of the effects of resonant waves, for which there is 
no ambiguity. It is however worth keeping in mind that the introduction of the non-resonant 
waves is likely to result in significantly higher maximum energies at the early stages of the 
Sedov- Taylor phase. 

Our complete calculations, including the non-linear dynamical reaction of the accelerated 
particles, the resonant amplification of ma gnetic field and the dy namical reaction of the field 



itself have been carried out as discussed by lCaprioli et al.l (J2008bl ). The results are illustrated 
in Figs. [2JH1 In the left panels we plot the maximum momentum (p max ), the shock velocity 
(V) and the two compression factors (R su b and Rtot) as functions of time. The right panels 
show the acceleration efficiency and the escape flux normalized to PoKl an d (l/2)poKi 
respectively, and the strength of the downstream magnetic field in units of 10 3 pG. The 
maximum momentum and the shock speed are in units of 10 6 m p c and 10 8 cms -1 respectively. 
The three figures refer to the following sets of parameters: no = 0.1cm -3 , Bq = lfiG (Fig. [2]), 
n = 0.1cm~ 3 , B = 5pG (Fig. ED and n = 0.03cm" 3 , B = lfiG (Fig. gj). 

The maximum momentum is determined at each time by requiring that the diffusion 
length in the upstream section equals 0.1R sh (t). The quantities p m ax{t)i R S ub{t) and Rtat(t) 
are all outputs of the non-linear calculations at the time t. The first point in time in all 
figures corresponds to the beginning of the Sedov phase. The time at which the Sedov- 
Taylor expansion begins was determined assuming E51 = 1 and that the mass of the ejecta 
is M e j = 5Mq. The other relevant parameters are the temperature of the ISM in which the 



10 4 K, and the momentum threshold for 



SNR is expanding, for which we assumed T 

particles to be injected into the accelerator, which was chosen to be p in j 

with ipi n j = 3.8 and T 2 the temperature downstream of the shock (see lCaprioli et al.l ( l2008bl ) 

for details). 

Some general comments are in order: one may notice that the total compression factors 



obtained in our calculations are always lower than ~ 10. This is unique 



Caprioli et al. 



y due to the dynam- 
fl2008ah the effect of 



ical reaction of the amplified magnetic field. As shown by 
the amplified field on the plasma compressibility is relevant whenever the magnetic pressure 
becomes comparable with the thermal pressure of the background plasma upstream. The 
consequent decrease in the compression ratios allows us to be consistent w ith the values that 
have been inferred from observations of a few SNRs (I Warren et al.l (120051 )). 

The highest momentum of accelerated particles, as expected, is reached at the beginning 
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Figure 2. Left panel: time dependence of the maximum momentum of accelerated particles (solid curve labeled as P m ax,6) in 
units of 10 6 m p c, of the shock velocity in units of 10 s cm/s (solid curve labeled as Vg), of the compression factor at the subshock 
R su t, (dashed curve) and of the total compression factor Rtot (dot-dashed curve), during the Sedov- Taylor phase. Right panel: 
time dependence of the magnetic field strength downstream of the subshock in units of 10 3 fiG (dot-dashed curve), of the 
escape flux normalized to poV^/2 (dashed curve) and of the cosmic ray pressure at the shock location normalized to poV\. 
The magnetic field strength, Bq, and the number density of the background plasma, no, at upstream infinity are taken to be 
Bo = 1 fJ'G and no = 0.1 cm -3 . 



of the Sedov- Taylor phase and is of order ~ 10 6 GeV (about the knee). This should be 
considered as a lower limit to the maximum energy reached at that time, since we have 
decided not to include the non-resonant channel of magnetic field amplification, which is 
very efficient when the shock velocity is large. During the following expansion, the time- 
dependence of Pmax is reasonably well approximated by p ma x oc £ -1 / 2 , in agreement with 
Eq. [MJ since £ c (£) is roughly constant (see solid curve on the right panels of Figs. EHUand 
discussion below). 

A crucial ingredient in calculating the maximum energy at a given time is the strength 
of the magnetic field. The magnetic field intensity in the downstream plasma is plotted in 
the right panels (dot-dashed line) for the three cases considered here. The typical values are 
between ~ a few— 10/iG at late times and ~ 30 — 100/iG at the beginning of the Sedov phase. 
After the first few thousand years, the dependence on time is not far from SB2 oc t _3//10 , as 
would result from Eq. [321 using the additional information on the approximate constancy 
of £ c (£) and R su b (dashed line in the left panels of Figs. I^JtoHJ). Despite these resemblances, 
the situation to which the plots refer is more complicated than that described by Eq. 1321 
where a number of effects have been neglected, first among these the presence of a precursor 
which evolves with time (Rtot is changing as can be seen from the dot-dashed curve on the 
left of Figs. [2] to H]) and the time varying adiabatic compression it entails. 

The right panels also show the acceleration efficiency (solid line) and normalized escape 
flux (dashed lines). One should notice that even when the acceleration efficiency is very high, 
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Figure 3. Same as Fig. [2]but for f?o = 5 fiG and no = 0.1 cm 
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of order ~ 50 — 60%, the escape energy flux never exceeds ~ 30%. As discussed above, this 
latter quantity should be the one that is more directly related to the cosmic ray energetics 
in the Galaxy, at least at the highest energies, while the former is more relevant for the 
generation of secondary radiation due to cosmic ray interactions in the acceleration region. 
The acceleration efficiency and the normalized escape flux initially increase with time 
during the Sedov- Taylor expansion phase. This behaviour is related to an analogous trend 
of the shock modification, as can be clearly seen from the time-dependence of R to t (dash- 
dotted curve in the left panel of Figs. EM]). In fact, at the beginning of the Sedov phase, 
the amplified magnetic field is at a maximum and its dynamical reaction on the shock is 
so strong that the acceleration efficiency is reduced. As soon as the magnetic field strength 
starts decreasing, the shock modification increases, and £ c and F esc with it. Notice, however, 
that this does not mean that the actual cosmic ray pressure and escape flux increase, because 
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£ c and F esc are normalized to PoV^ h (t) and p V^ h (t)/2 respectively, and both decrease with 
time rather quickly. 

At later times, both £ c and F esc start decreasing, with the latter showing a more rapid 
decline than the former. This is due to the fact that the shock is slowing down and progres- 
sively becoming unmodified: the maximum momentum is decreasing and the spectrum of 
accelerated particles is steepening. Recalling again that the plots show normalized quanti- 
ties, one gathers that the decline with time of cosmic ray pressure and escaping energy flux 
is quite dramatic in this phase. 

6 CONCLUSIONS 

The escape of accelerated particles from the shock region towards upstream appears in basi- 
cally all approaches to non-linear particle acceleration at shock fronts. In kinetic approaches 
that solve the relevant equations assuming quasi- st at ionarity at any given time the escape 
flux appears either as a corrective term in the equation of energy conservation between 
downstream and upstream infinity, or as a result of imposing as a boundary condition that 
the distribution function of accelerated particles vanishes at a finite distance xq upstream 
of the shock. In this second case the escape flux is described by the non-vanishing spatial 
derivative of the distribution function at x . Such a quasi-stationary description of the ac- 
celeration process is meaningful only whenever the shock velocity decreases in time and the 
magnetization at the shock, as due to streaming of cosmic rays, also decreases in time. These 
two conditions lead to escape of the particles with momentum close to the highest achievable 
at any given time. This would be the case during the Sedov- Taylor phase of the expansion of 
a supernova shock. These approaches are therefore unfit to describe the phase of free expan- 
sion, where the maximum momentum increases with time and no escape towards upstream 
takes place. 

Fully time-dependent, numerical calculations apply to arbitrary situations and should 
be able to describe the increase of momentum during the free expansion as well as the es- 
cape of particles during the Sedov phase. However, in most cases these calculations do not 
include the time dependence of the magnetic field at the shock. This leads to an increases 
of the maximum momentum even during the Sedov phase. Probably for this reason there 
is no detailed discussion in the literature of the escape flux from a SNR in the context of 
time dependent calculations. The appearance of the escape flux in numerical time depen- 
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dent approaches would depend on the boundary conditions adopted in the specific case. If 
the distribution function is assumed to vanish at upstream infinity, the escaping particles 
should appear as a concentration of particles with large momenta (close to p m ax{t)) at large 
distances from the shock: since the acceleration box extends to infinity, there is no real es- 
cape. On the other hand, if the boundary condition f(xo) = is imposed at a finite distance 
from the shock, the escape flux (as a function of time) should appear automatically as a 
consequence of a non- vanishing Ddf(xo)/dx. It would be interesting to have confirmation of 
these expectations by running time dependent calculations with decreasing magnetization. 
The existence of th e escap e flux is not simply a mathematical nuisance: as discussed by 



Ptuskin fc Zirakashvilil (120051 ) it is the very reason why supernova remnants can potentially 



be the sources of cosmic rays up to the knee. If the particles we observe as cosmic rays were 
the ones advected towards downstream of the shock, the adiabatic losses suffered during the 
expansion of the shell would drive these particles towards lower en ergies. It is therefore of 



Ptuskin fc Zirakashvili 



the hi ghest importance to go beyond the test particle calculation of 

(120051 ) and achieve an understanding of the escape flux in the context of the fully non-linear 

theory of diffusive shock acceleration. 

In this paper we made a first step in this direction, by calculating the temporal evolution 
of the escape flux during the Sedov- Taylor phase of the expansion of a supernova in a 
medium with constant density. Our calculation included a recipe to describe the cosmic 
ray induced magnetic field, the dynamical reaction of the accelerated particles and the field 
itself. In addition to this we also discussed in some detail the limitations of the conservation 
equations as they are usually written, in that they do not allow to describe the magnetic 
field perturbations when they are not Alfven waves. 

We showed that the escape flux may involve between few and 10-30 % of the shock ram 
pressure, while the particle acceleration efficiency at the same time reaches 50-60 %. This 
means that the energetics inferred from observations of secondary radiation from a remnant, 
for instance in the form of gamma rays, may not be the relevant one for the origin of cosmic 
rays, since it is not related in a trivial way to the energy flux escaping the accelerator. 

The maximum energy up to which particles may get accelerated is reached at the begin- 
ning of the Sedov phase and is of order 10 15 eV if only resonant amplification of the field 
is included. E max might be somewhat larger for some SNRs that at the very beginni ng of 
the S edov- Taylor phase may experience the effect of non-resonant streaming instability (JBell 



20041 ). This mechanism provides extremely efficient field amplification as long as the shock 
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velocity is high, and hence is expected to play a very important role in the early Sedov- Taylor 
phase (and possibly during the free expansion phase). At present, non-resonant modes can- 
not be formally accounted for in the conservation equations, since a detailed description 
of the energy transfer between particles and waves is not available yet. We carried out all 
calculations in the simpler case of Alfven waves interacting with particles in a resonant way. 
The amplification of the magnetic field, by either resonant or non-resonant streaming 
instability, has profound implications on the escape flux of particles towards upstream of 
a shock, and therefore on the spectrum of cosmic rays we observe at the Earth. The most 
obvious consequence of the magnetic field amplification is that of allowing for higher values 
of the maximum energy of accelerated particles, as shown by our Eqs. [34] and [351 However 
large magnetic fields exert a dynamical reaction on the plasma leading to a reduction of 



the compression in the precursor. 



the pressure of the background gas (ICaprioli et al 



his happens whenev er the magnetic pressure exceeds 



spectrum of accelerated particles (ICaprioli et al.l l2008bl ) is reduced and at the same time 



2008a]). As a result, the concavity of the 



the escape flux at p ~ p ma x decreases. It follows that larger field implies larger p max but not 
necessarily larger escape flux, as shown in Figs. H] [2] and [3] (see the behaviour of the curves 
at early times). 

The escape of accelerated particles from a cosmic ray modified shock has profound im- 
plications for the origin of cosmic rays, which will be discussed in detail in a forthcoming 
paper. Here we want to emphasize some general points: 

1) the escape from upstream is the natural solution to the well known problem of ex- 
plaining how the highest energy particles (say, close to the knee energy) could escape the 
system without suffering substantial adiabatic energy losses; 

2) the magnetic field amplification is expected to switch from mainly non-resonant to 
mainly resonant at the beginning of the Sedov- Taylor phase. It can be easily understood 
that this may lead to peculiar changes in the spectrum of cosmic rays detected at the Earth, 
reflecting this transition; 

3) the flux of escaping particles, once integrated in time during the SNR evolution may be 
very different from the concave instantaneous spectrum which can potentially be observed 
in a SNR, for instance by looking at its multifrequency emission. This point is certainly 
relevant for the purpose of addressing the commonly raised point of how a concave spectrum 
of accelerated particles can reflect in an almost perfect power law over many orders of 
magnitude; 
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4) there is a further complication of all the picture, due to the acceleration of nuclei at 
energies that may be expected to scale as the charge of the nucleus (in the case of Bohm 
diffusion). Any calculation of the flux of single chemical species observed at the Earth must 
take these complex effects into account. 
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